Measurement of charmless B Decays to r\K* and r\p 
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Yonsei University, Seoul 

We report measurements of branching fractions and CP asymmetries for B —> r\K* and B — > np 
decays. These results are obtained from a 414 fb -1 data sample collected at the T(4S) resonance 
with the Belle detector at the KEKB asymmetric-energy e + e~ collider. We measure the following 
branching fractions: B{B° -» r)K*°) = (15.2 ± 1.2 ± 1.0) x 1(T 6 and B{B + -> T)K* + ) = {19.2>±\% ± 
1.5) x 10 -6 , where the first error is statistical and the second systematic. We also find a 2.7a 
excess in the B + -> np + mode and measure B{B + -> rjp + ) = (4.1±^;| ± 0.4) x 10~ B < 6.5 x 10~ 6 
at 90% confidence level. For B° — > np decays, we determine the upper limit B(B a — > ??p ) < 
1.9 x 10" 6 at 90% confidence level. The partial rate asymmetries are Acp\r}K*°) = 0.17±0.08±0.01, 
A C p{nK*+) = 0.03 ± 0.10 ± 0.01, and Acp(vP + ) = — 0.04±g;|| ± 0.01. 

PACS numbers: 13.25.Hw,14.40.Nd 



INTRODUCTION 

Charmless hadronic B decays play an important role 
in understanding CP violation in the B meson system. 
The decays B — » t]K* and B — > r\p are key examples. In 
the standard model (SM), penguin (tree) diagrams are 
expected to dominate in B — > r\K* (B — > r\p) decays 
(Fig. [J). The large branching fraction for B — > r)K* 
compared to that for B — > r\K P, 0, Q can be explained 
qualitatively in terms of the interference between non- 
strange and strange components of the r\ meson, but is 
higher than recent theoretical predictions 0, H, [a, 0] • In 
a similar vein, the larger measured branching fraction 
for charged (B + — * r]K* + ) versus neutral (B° — * r]K*°) 
decays may suggest an additional SU(3)-singlet contribu- 
tion IH or constructive interference between SM pen- 
guin and tree amplitudes or between SM and new physics 
penguin amplitudes. Throughout this paper, the inclu- 
sion of charge-conjugate modes is implied unless stated 
otherwise. 

In the standard model, direct CP violation (DCPV) 
occurs in decays due to interference between two (or 
more) amplitudes that have different strong and weak 
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FIG. 1: Feynman diagrams for B + — > rjK* + and np + decays. 
The corresponding neutral decays are similar except that the 
spectator quark becomes a d and (b) and (c) diagrams do not 
exist. 



phases. The partial rate asymmetry can be written as 

CPK ' r(s —* f) + r(5 — > /) 

_ 2|Mi||M 2 | sin A(5sin A0 

~ |Mx| 2 + \M 2 \ 2 + 2|Mi||M 2 |cosA5cosA(/>' 



(1) 
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where B and / are the CP-conjugate states, and AS (A</>) 
is the difference of the strong (weak) phases between am- 
plitudes Mi and M 2 . Here, the amplitude M 1 (M 2 ) rep- 
resents the sum of the amplitudes from penguin (tree) 
diagrams having a common weak phase. The asymme- 
try will be sizable when the two type of amplitudes are 
of comparable strength with significant phase differences. 
However, B — ► rjK * and B — > rjp decay rates are expected 
to be dominated by penguin diagrams for rjK* and tree 
diagrams for 77/9, so Acp is expected to be small. On the 
other hand, amplitudes arising from new physics may in- 
terfere with these SM amplitudes to generate a sizable 
Acp value. The experimental results [3| suggest that 
DCPV is small, albeit with large statistical uncertain- 
ties. 



DATA SET AND APPARATUS 

This analysis is based on a data sample collected at 
the T(4S) resonance with the Belle detector at the 
KEKB [10| accelerator. The data sample corresponds to 
an integrated luminosity of 414 fb -1 and contains 449 x 
10 6 BB pairs. 

The Belle detector is designed to measure charged 
particles and photons with high efficiency and preci- 
sion. Charged particle tracking is provided by a sil- 
icon vertex detector (SVD) and a central drift cham- 
ber (CDC) that surround the interaction region. The 
charged particle acceptance covers the laboratory polar 
angle between = 17° and 150°, measured from the z 
axis that is aligned anti-parallel to the positron beam. 
Charged hadrons are distinguished by combining the re- 
sponses from an array of silica aerogel Cherenkov coun- 
ters (ACC), a barrel- like array of 128 time-of-fhght scin- 
tillation counters (TOF), and dE/dx measurements in 
the CDC. The combined response provides K/n sepa- 
ration of at least 2.5a for laboratory momentum up to 
3.5 GeV/c. Electromagnetic showers are detected in an 
array of 8736 CsI(Tl) crystals (ECL) located inside the 
magnetic volume, which covers the same solid angle as 
the charged particle tracking system. The 1.5-T mag- 
netic field is contained via a flux return that consists of 
4.7 cm thick steel plates, interleaved with resistive plate 
counters used for tracking muons. Two inner detector 
configurations were used. A 2.0 cm beampipe and a 3- 
layer silicon vertex detector were used for the first sample 
of 152 x 10 6 BB pairs, while a 1.5 cm beampipe, a 4-layer 
silicon detector and a small-cell inner drift chamber were 
used to record the remaining 297 x 10 6 BB pairs [ill ]. 

We calculate the acceptance and study backgrounds 
using Monte Carlo (MC) simulation. For these simula- 
tion studies, the signal events, generic b — > c decays and 
charmless rare B decays are generated with the EVT- 
GEN |12j event generator. The continuum MC events 
are generated with the e + e~ — > 7* — > qq process in the 



JETSET [13J generator. The GEANT3 [lj package is 
used for detector simulation. 



EVENT SELECTION AND RECONSTRUCTION 

Hadronic events are selected based on the charged 
track multiplicity and total visible energy sum, which 
give an efficiency greater than 99% for BB events. All 
primary charged tracks are required to be consistent with 
coming from the run-dependent interaction point within 
±2 cm along the z axis and within ±1.5 cm in the trans- 
verse plane. Particle identification (PID) is based on 
the likelihoods Ck and £ T for charged kaons and pi- 
ons, respectively. These likelihoods are calculated from 
CDC, TOF, and ACC information. A higher value of 
£k/(£tt + £k) indicates a more kaon-like particle. In 
this analysis, PID cuts are applied to all charged particles 
except those associated with Kg — > ir + ir~ decays. Un- 
less explicitly specified, the PID cuts are £k/(£tt + £k) 
> 0.6 for kaons and < 0.4 for pions. The corresponding 
efficiencies are 85% for kaons and 89% for pions; 8% of 
pions are misidentified as kaons and 11% of kaons are 
misidentified as pions. 

We form ir° candidates from photon pairs with an 
invariant mass between 118 MeV/c 2 and 150MeV/c 2 
(2.5cr). The photon energies must exceed 50 MeV, and 
the 7T° momentum in the center-of-mass (CM) frame 
must exceed 0.35 GeV/c. Kg candidates are recon- 
structed from pairs of oppositely charged tracks whose 
invariant mass lies within ±10MeV/c 2 (2.5a) of the Kg 
meson mass. We also require that the vertex of the Kg 
be well-reconstructed and displaced from the interaction 
point, and that the Kg momentum direction be consis- 
tent with the Kg flight direction. 



77 Meson Reconstruction 

Candidate 77 mesons are reconstructed in the 77 — > 77 
and 77 — * 7r + 7T~7r modes. If one of the photons from 
the former 77 decay mode can be paired with another 
photon and have a reconstructed 77 mass within 3a of 
the 7T° meson mass, the 77 candidate is discarded. We 
relax the PID requirement for charged pions from the 
latter 77 decay mode to £k/(£tt + Lk) < 0.9. Can- 
didate 77 mesons are required to satisfy the following 
mass selections: 500 MeV/c 2 < M 77 < 575 MeV/c 2 and 
537 MeV/c 2 < M n+n -„o < 557 MeV/c 2 , where the recon- 
structed mass resolutions are 12 MeV/c 2 for 77 — > 77 and 
3.5 MeV/c 2 for 77 — ► 7r + 7r~7r°. When reconstructing the 
B meson candidate, the momentum of the 77 candidate 
is recalculated by applying the 77 mass constraint. The 
77 — > 77 candidates must satisfy | cos6>*| < 0.90, where 8* 
is the angle between the photon direction and the direc- 
tion of the CM frame in the 77 rest frame; this requirement 
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suppresses soft photon combinatorial and B — > K*j feed- 
across backgrounds. 



K* and p Meson Reconstruction 

K*° candidates are reconstructed from K + ir~ and 
Kgir pairs, while the K* + candidates are reconstructed 
from K + ir° and Kgir + pairs. These candidates are re- 
quired to have reconstructed masses within ±75 MeV/c 2 
of the nominal value [HI]. Candidate p° (p + ) mesons are 
reconstructed from ir~n + (7T°7r + ) pairs. Each combina- 
tion is required to have a reconstructed mass within ±150 
MeV/c 2 of the nominal value [l5j |. 



B Meson Reconstruction 

The B meson candidates are reconstructed from r/K*°, 
rjK* + , r/p°, and r/p + combinations. They are char- 
acterized by the bcam-energy-constrained mass M^ c = 
^/_E 2 cam /c 4 — \Pb/c\ 2 and the energy difference AE = 
E B - E hcam , where E heam = 5.29 GeV, and P B and E B 
are the momentum and energy, respectively, of the B 
candidate in the CM frame. We define the fit region in 
the M hc -AE plane as M hc > 5.2GeV/c 2 and \AE\ < 
0.25 GeV. We define the signal region as the overlap of 
the bands M bc > 5.27GeV/c 2 and \ AE\ < 0.05 GeV. 

From signal MC, we find that 8 — 10% of the events 
contain multiple B candidates. Only one B candidate 
per event is retained for the likelihood fit. If there are 
multiple r\ candidates, we choose the one with the small- 
est x 2 °f the fit with a mass (vertex and mass) constraint 
to the kinematics of the r\ meson in the case of rj — > 77 
(rj — > 7T + 7r~7r ) decays. Among the B candidates made of 
the same rj candidate, we choose the one with the smallest 
vertex x 2 m the cases of K*° — > K + tt~ or p° — > 7r + 7r~; 
or then the one with the mass closest to nominal in the 
case of K* + — > K°ir + ; or then the tt° mass closest to the 
nominal in all other cases. 



BACKGROUND SUPPRESSION 

The dominant background for exclusive two-body B 
decays comes from the e + e~ — > 7* — > qq continuum 
(<7 = u, d, s, c), which has a jet-like event topology in 
contrast to more spherical BB events. The other ma- 
jor backgrounds involve feed-across from these and other 
charmless B decays. The background from b — ► c decays 
has a small impact because the M^ c and AE distributions 
do not peak in the signal region. In this analysis, the 
fit does not distinguish non-resonant B — > rjK tt decays 
from B — > r/K* decays, since they have the same M^ c 
and AE distributions. The non-resonant contribution is 



estimated and subtracted independently using the Kir 
invariant mass distributions of the fitted -B-decay yields. 

Continuum Background 

Signal and continuum events are distinguished in two 
steps. Here, all the variables are calculated in the CM 
frame. First, we require |cos6>t| < 0.9, where 6t is de- 
fined as the angle between the 77 direction of a B candi- 
date and the thrust axis from all particles in the event 
not associated with that B candidate. This retains 90% 
of signal and removes ^56% of continuum. Second, a 
likelihood C s (£ c ) for signal (continuum) is formed from 
two independent variables — cos0 B , where 9 B is the po- 
lar angle of the B candidate momentum direction, and a 
Fisher discriminant [l6| T — a ■ R that combines seven 
event shape variables: cos 9t, Si (the sum of the magni- 
tudes of the momenta transverse to the rj direction for all 
particles more than 45° away from the 77 axis, divided by 
the sum of the magnitudes of the momenta of all particles 
not from the candidate B meson [13]): and the five mod- 
ified Fox- Wolfram moments [HI Rf, R° 2 °, R%°, and 
R1° . The Fisher discriminant's weight vector a is deter- 
mined to maximize the separation between signal events 
and continuum background using MC data; these Fox- 
Wolfram moments are used since they are not correlated 
with Mb c . The likelihood ratio 1Z = C s / (C s +C C ), which 
peaks near one for signal and near zero for continuum, is 
used to distinguish signal from continuum. 

The distribution of 1Z is found to depend somewhat 
on the event's B flavor tagging quality parameter r [l9| . 
which ranges from zero for no flavor identification to 
unity for unambiguous flavor assignment. We partition 
the data into three r regions, r < 0.5, 0.5 < r < 0.75, and 
r > 0.75. In each r region, the optimal cut on 1Z is de- 
termined by maximizing the significance Ns /y/Ng + N B , 
where Ns and N B are the retained number of signal and 
continuum background events selected in MC samples. 
For cut optimization studies, we assume the branching 
fractions of 2.0 x 10~ 5 for rjK*, 5.0 x 10" 6 for r)p+, and 
1.0 x 10~ 6 for rjp°. For B° -> rjK* ^ -> 77, iT* -> 
K + n~ decays, a typical cut of 1Z > 0.4 is ~ 87% ef- 
ficient for signal and removes ~ 67% of the continuum 
background for data in the region r < 0.5, and a cut 
of 1Z > 0.2 is ~ 94% efficient for signal and removes 
~ 53% of the continuum background for data in the re- 
gion r > 0.75. 

Backgrounds from B Decays 

B — * K*(p)j is the dominant charmless B-decay back- 
ground for B — * r]K*(p), rj — > 77 decays. The 77 — * 77 
selection, |cos6**| < 0.90, removes 85% of this back- 
ground. To further suppress it, we pair each photon 
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from the rj — > 77 candidate with the K* or p candi- 
date and reject those events where M^ c > 5.27GeV/c 2 
and -0.2 GeV < AP < 0.1 GeV. We thus remove 96% of 
this background and retain 93% of the signal events. For 
P — ► rjp decays, a measurable contribution from other 
charmless P and b — > c decays remains (see Table IJ). 
The contributions of these backgrounds are taken into 
account in the analysis. 

TABLE I: Estimated B — decay background contributions in 
the fit region to B — > r\p from 6 — > c (iVbc), charmless P 
decay (AT r ), and jyTf* feed-across (A^ ee d) and measured yields 
from all sources (N) but dominated by residual continuum 
background, after application of the cos#t and 1Z cuts. Nbc 
and N T are estimated from MC samples. 



Mode 


N bc 




N 


V-f~tP° 


62 


81 


17 


2931 


Vttitit® P 


67 


27 


5 


1063 


V-nP + 


148 


74 


3 


4169 


VnnnO P + 


76 


22 


1 


1809 



ANALYSIS PROCEDURE 

Signal yields are obtained using an extended unbinned 
maximum likelihood fit to the Mbc and AP distributions 
(2-D ML) for events that satisfy the cos 8t and 1Z re- 
quirements. 

For N input candidates, the likelihood is defined as 

(jVs+A^+Wbc + Nr+Wfccd) JL 

L = W[ Y[(N S P Si + 

»=1 

N m p w, + ^bc^bc, + N r P n + iV fccd P fccdi ), (2) 

where Ps^ , P qq . , Ph Ci , and Pf e e di are the probability 
density functions for event i, with measured values Mhc,i 
and AEi, to arise from signal, continuum background, 
b — > c background, charmless B decay background, and 
feed-across background, respectively. The small yields 
iVbc, N T) and ATf eed are fixed from the MC analysis. 

The continuum, b — > c and charmless P-decay back- 
ground AP probability density functions (PDF) are mod- 
eled by second- or third-order polynomial functions. The 
continuum and b — + c background components in Mb c 
are modeled by a smooth function [20|. To account 
for the peaking behavior of Mt, c in the signal region 
from charmless P decay backgrounds, we use the sum 
of two bifurcated-Gaussian functions to model the dis- 
tributions. The bifurcated Gaussian combines the left 
half of a wide-resolution Gaussian with the right half 
of a narrow-resolution Gaussian, both having a common 



mean. For P — > r\p decays, the M^ c and AP distribu- 
tions from r]K* feed-across will behave like signal with a 
AP shift of -50 MeV. The PDF shape for each contribu- 
tion is determined from MC. The first-order coefficient of 
the continuum-background AP polynomial and the pa- 
rameters of the Mb c function are allowed to float in each 
fit. 

For the signal AP distribution, we combine two 
bifurcated-Gaussian functions. The first accounts for 60- 
80% of the total area and the wider second models the 
low-energy tail. M^ c is weakly correlated with AP, so 
we construct separate bifurcated Gaussians for M^ c in 
the three ranges |AP| < 0.05 GeV, 0.05 GeV < |AP| < 
0.1 GeV, and 0.1 GeV < \ AE\ < 0.25 GeV. The parame- 
ters of these functions are estimated from MC first, then 
calibrated with a large control sample of P + — > D°ir + , 
U° -> K+tt-tt decays. 

For decays with more than one sub-decay process, the 
final results are obtained by fitting the sub-decay modes 
simultaneously with the expected efficiencies folded in 
and with the branching fraction as the common output. 
The statistical significance (E) of the signal is defined as 
y— 2 ln(Lo/£ m ax), where Lq and L max denote the likeli- 
hood values for zero signal events and the best fit num- 
bers, respectively. 

The 90% confidence level (C.L.) upper limit xgo on the 
signal yield is calculated from the equation 

/? 90 L(x) dx 

J oo r , =90%- 

Jo L(x)dx 

To incorporate the systematic uncertainty in the calcu- 
lation of £90, the likelihood function is smeared with a 
Gaussian function with the resolution from the system- 
atic uncertainty. That smeared likelihood function is also 
used to calculate the significance of the signal including 
the systematic uncertainty. 

MEASUREMENTS OF BRANCHING 
FRACTIONS 

Efficiencies and Corrections 

The overall reconstruction efficiency e is first obtained 
using MC samples and then multiplied by PID efficiency 
corrections obtained from data. The PID efficiency cor- 
rection is determined using D* + — > D°ir + , D° — > K~ir + 
data samples. Other MC efficiency corrections are deter- 
mined by comparing data and MC predictions for other 
well-known processes. The charged-particle tracking effi- 
ciency correction is studied using a high-momentum 77 
sample and is determined by comparing the ratios of 
rj — > 7r + 7T~7r° to rj — > 77 in data and MC. The same 
high-momentum rj sample is also used for ir° recon- 
struction efficiency corrections by comparing the ratio 



of 77 — > urn it to 77 — > 77 between the data and MC 
sample. The Kg reconstruction efficiency is verified by 
comparing four K*(892) decay channels (K + ir~, K + ir°, 
K s tt + , K s w°) in inclusive K* and exclusive B — > J/ijjK* 
samples. The 1Z cut efficiency correction is determined 
using B + — ► D°ir + decays. For 77 and K* reconstruction 
and mass cuts, we use the high-momentum 77 and K* 
sample for the efficiency correction studies. The above 
studies show good agreement between data and MC; the 
reconstruction and selection efficiencies differ by about 
2%. The P1D, 7T°, 77 and K* reconstruction efficiency 
corrections are applied and the systematic uncertainties 
are also obtained from the above studies. 



TABLE II: Summary of results for each channel listed in the 
first column. The measured signal yield (Ns), reconstruction 
efficiency (e), total efficiency (etot) including the secondary 
branching fraction, statistical significance (E) and measured 
branching fractions are shown. Uncertainties shown in second 
and sixth columns are statistical only. For the final combined 
branching fractions, corrections for contributions from non- 
resonant or higher resonance components have been applied. 
The total systematic uncertainties are given, and the com- 
bined significances include the systematic uncertainties. 



Mode 


N S 


e(%) 


etot(%) 


E 


£>(10 -6 ) 


'Innn ' v k+tv~ 


336.2 + ^'2 
93.4±i*:i 
20.1+1? 


16.9 
9.8 
2.1 
1.3 


4.4 
1.5 
0.27 
0.098 


14.2 
8.7 
3.6 
2.6 


16.9 ± 1.5 

14.13:? 

16.7± 6 5 .6 
21.6tft B 


riK*° 








15.7 


15.2 ± 1.2 ± 1.0 


Vl-y K-k+tt 

Vititit K f(0„+ 


24.i±?:S 

120.3t]|;2 

00 o+ 7 - 3 


6.7 
4.2 
4.5 
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0.88 
0.32 
1.2 
0.38 


6.1 
3.5 
10.1 
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16.9^i 
22.6+29 
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nap 

^?7T7r7r^ P 
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8.9 
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1.2 


2.1 
0.2 


1 nc+0.73 

1 - ZO -0.67 
U.-W _o.66 


VP° 








1.3 


0.84lol5? ±0.19 


ni~,p + 

^?7T7r7r^ P 


38.lll«; 2 

i5.8±S:g 


5.5 
3.50 


2.2 
0.79 
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4.1±i;| ±0.4 
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FIG. 2: Projections on Mb c (for the signal slice in AE) and 
AE (for the signal slice in Mbc) for r)K*° (a,b) and r/K* + (c,d) 
with the expected signal and background curves overlaid. The 
shaded area represents 77 — > tt + iy~tt° decays. 




5.25 

M K (GeV/c 2 ) 
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FIG. 3: Projections on Mf, c (for the signal slice in AE) and 
AE (for the signal slice in Mbc) from 2-D ML fit results for rjp° 
(a,b) and rjp + (c,d) with the expected signal and background 



curves overlaid, 
decays. 



The shaded area represents 77 — » tt 7r 



Figure [5] shows the projections of the data and the fits 
onto Mbc (for events in AE signal slice) and AE (for 
events in the Mt> c signal slice) for the B — ► 77^* decays, 
while Fig. [3] shows the corresponding projections for the 
B —>■ rjp decays. 



Non K* (892) Components 



Fit Results 

The fitted signal yields and branching fractions are 
shown in Table [TlJ Several consistency checks are made, 
including tighter TZ cuts as well as 1-D ML Mbc and AE 
fits, and they are all found to be consistent. The total 
observed yields are N, lK ,o = 459.2+j^ for B° -> r/K*°, 
N nK ,+ = 253.4+J5-5 for B+ -► r)K*+, iV„„o = 20.4 



12.2 
-11.0 



for B 







VP 



and 2V_ 



53.9ti 8 7 .t 



for 5 



The background-subtracted fsT* helicity distributions 
within the Mb c and AE signal regions (Fig. [4]) are con- 
sistent with the expectation from rjK* decays, indicating 
no significant S-wave or higher resonance contribution in 
the K* mass region. The K*°(K *+) helicity angle (9 hcl ) 
is the angle between the vT _ (7r°, K°) direction and the 
opposite of the B direction in the K* rest frame. The 
binned x 2 P er degree of freedom is x 1 /N = 9.7/10 for 
K*° and x 2 /N = 3.6/10 for K*+. 

We use a 2D-ML fit to the Kir invariant mass distri- 
butions to evaluate the small S-wave or higher resonance 
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cays. These corrections are applied to the final branching 
fraction measurements of B — > r\K * . 



For B^ 



r\p decays, we examine the properties of 



the p + candidates through 2-D ML fits in bins of iy + tt° 
invariant mass and p + hclicity. Although statistically 
limited, a p + mass peak and a polarized p + helicity dis- 
tribution are observed (Fig. [6]) and are consistent with 
the expectation from rjp + . Due to the limited statis- 
tics for B + — ► r)p + decays, a larger systematic error for 
the non-resonant or higher resonance contributions is as- 
signed with no corrections applied. 



FIG. 4: Distributions of the K* helicity for the (a) K*° and 
(b) K* + modes in case of B — * rjK* decays. The overlaid 
histograms represent the distributions from MC normalized 
by the 2-D fit results. 




1.27 1.46 1.65 0.70 0.8 
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FIG. 5: Fitted yields vs. the Kn invariant mass for the (a) 
(Kit) and (b) (K-k) + modes. The overlaid functions are the 
results of a binned \ 2 fit- The dashed line represents the con- 
tribution from the D-wave ^1(1430), and the dot-dash line 
represents the LASS S-wave parameterization. LASS parame- 
terization parameters, widths of P-wave and D-wave functions 
are allowed to float in the fits. 



SYSTEMATIC ERROR 

Systematic errors, enumerated in Table HlH arise from 
efficiency corrections, non-resonant corrections and fit- 
ting. The main sources of uncertainties in the efficiency 
corrections are from the reconstruction of low- momentum 
charged tracks, low-energy photon finding, and the 1Z cut 
efficiency. The systematic errors include contributions of 
1% for 1Z cuts, 1% per reconstructed charged particle, 
0.5% for each charged particle identification, 4% for ir° 
reconstruction, 4.5% for Kg reconstruction, and 2% for 
r\ reconstruction with r\ — > 77. We use B + — ► D°ir + 
decays to estimate the uncertainties in the signal PDF's 
used for fitting M^ c and AE by comparing the mean 
and the width of the M^ c and AE distributions be- 
tween the B + — > D°ir + data and the MC sample. In 
Table IIII1 "Fit" means the systematic uncertainty from 
the PDF function modeling. U B S " means the systematic 
uncertainty from the braching fractions of i] and K*(p) 
decays, which is obtained from the PDG tables [la ]. 
"Non-resonant" means the systematic uncertainty from 
the non-resonant or higher resonant contributions in the 
K*(p) mass window regions, which is obtained from the 
studies of a 2D-ML fit to the Kir invariant mass dis- 



Kir contaminations in the K* mass region. A clear ex- 
cess in the higher Kit invariant mass region is observed 
(Fig. [5]) . To estimate these contributions, we fit the dis- 
tributions with a P-wave relativistic Breit-Wigner func- 
tion for the if* (892), a D-wave relativistic Breit-Wigner 
function for the K% (1430) resonance and an ad hoc func- 
tion for the S-wave contribution. Several functions are 
used to model the S-wave contribution in the K ir mass 
region 1.0 GeV/c 2 to 1.5 GeV/c 2 , including a if,* (1350) 
resonance [2 lj. the LASS distribution 2l|, and a threshold 
function. The LASS distribution contains a non-resonant 
S-wave background function interfering with an S-wave 
K q (1430) resonance. Figure [5] shows one example of our 
fitting results with the S-wave contribution modeled by 
the LASS distribution. Based on these studies with vari- 
ous S-wave functions, the non-resonant Kir contributions 
are (5.6 ±3.0)% for r]K* Q and (5.0 ±3.0)% for r]K*+ de- 



(a)rip + 



_i 1 i_ 




M„„ (GeV/c 2 ) 



FIG. 6: Fitted yields vs. (a) p + helicity and (b)7r + 7r° invariant 
mass from B — > rjp + decays. The overlaid histograms are 
expected distributions from MC and normalized by the 2-D 
fit results with j^/jV = 6.8/8 for (a) and X*/N = 13.6/11 for 
(b). 
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tribution. For MC estimated b — > c and charmless B 
decay backgrounds ('Wb c " and "ZV r ", respectively), we 
vary the estimated yields by ±50% and refit the data. 
The difference between the resulting signal yield and the 
nominal value is taken as an additional systematic error. 
The overall relative systematic errors are 6.5% for r/K* , 
7.5% for r]K*+, 22.1% for r) P ° and 9.6% for r?p+. 



TABLE III: Relative systematic errors for r)K* and r/p. The 
unit is in percent (%). 



Contribution 




riK* + 


VP° 


VP + 


charged track./ Kg reconstruction 


2.9 


4.4 


3.3 


1.8 


7r°/?7 — ► 77 selection 


2.8 


3.5 


3.2 


4.8 


?7 mass window 


2.0 


2.0 


2.0 


2.0 


K*(p) mass window 


2.0 


2.0 


2.0 


2.0 


PID correction 


1.3 


1.0 


1.6 


0.8 


1Z requirement 


1.0 


1.0 


1.0 


1.0 


Fit 


1.7 


1.8 


7.5 


4.2 


iVbc, N T 


1.2 


1.2 


7.5 


0.4 


r}K* feed-across 






17.4 


1.1 




1 


1 


1 


1 


B s 


1.0 


1.0 


1.4 


0.8 


Non-resonant 


3.0 


3.0 


6.0 


6.0 


Total 


6.5 


7.5 


22.1 


9.6 



Acp MEASUREMENTS 

We measure Acp for B — > r\K* and B + — > rjp + . 
To account for the wrong-tag fraction w, the true value 
of Acp is related to the measured A^p via A'qp = 
(1 — 2w)Acp- Among the decay modes we study, only 
those in which the Ac p values are determined by low mo- 
mentum charged pions have a significant w: the wrong- 
tag fractions for K*+ -► K°tt+ is - 1.5% for r]K* + 
decays and ~ 2.0% for rjp + decays, while other de- 
cays have w < 0.1%. Since the result for r)K* + is 
obtained from a simultaneous fit to all four sub-decay 
modes with roughly equal statistics for K* + — > K + w 
and K* + — > K°ir + , the wrong-tag effect for Acp is 
less than 0.7%. Therefore, the only mode where we 
apply a correction due to the wrong tag fraction is 
B + — > r)p + , where we estimate w — 2%. To incorpo- 
rate the CP asymmetry in the fit, the coefficients of the 
signal and continuum background PDF's in the likeli- 
hood are modified as follows: Ns — > |A^s(l — lA^p) 
and N qq -> ^N gq (l - qA C p, qq ), where q = +1(-1) for a 
B(B) meson tag and A c °p,Acp,qq are the Acp outputs 
for signal and continuum, respectively. The results are 
A°£ P {j}K*°) = 0.17 ±0.08, A° c r >p(r 1 K*+) = 0.03 ± 0.10 
and A&(r,p+) = -O.Mjftig. 

Since the systematic errors in the reconstruction of 



the r) candidates and the number of BB events cancel 
in the ratio, the systematic uncertainty on Acv comes 
mainly from the charge asymmetry in the identification 
of charged kaons and the fitting PDF's. To estimated the 
fitting-PDF systematic uncertainty, we apply the same 
procedures as in the branching fraction measurements. 
The relative systematic errors from fitting PDF's are es- 
timated to be 3% for yyAT* , 13% for i]K* + , and 27% for 
i]p + . The efficiency asymmetry for the PID of charged 
kaons is 0.01 in absolute value. 

SUMMARY 

In summary, we report measurements of the exclusive 
two-body charmless hadronic B — > i]K* and B — > r\p 
decays with high statistics. Our results are consistent 
with previous measurements [H, 0| and confirm that the 
branching fractions for B° — * r)K*° and B + — > r\K* + 
are large. The branching fractions obtained are B(B° — > 
v K* a ) = (15. 2± 1.2± 1.0) x 10" 6 , and B(B+ -> r)K*+) = 
(19. 3+^g ± 1.5) x 10" 6 , where the first error is statisti- 
cal and the second systematic. Our measurements in- 
dicate that the branching fraction for B + — > r\K* + is 
1.4(7 higher than that for B° — * r/K* . A 2.7cr excess 
is seen for B + — * rjp + decays. The branching frac- 
tion and 90% C.L. upper limits for B — > r/p decays are 
B(B+ -> rip+)= (4.ll^±0.4) x 10- 6 (< 6.5 x lO" 6 ) and 
B(B° 7?p°)< 1.9 x 10~ 6 . The measurements of the 
B — > r/p branching fractions are consistent with theoret- 
ical predictions jj 0,0,0 ■ 

We have measured the direct CP asymmetry in the 
B — > rjK* and B + — * r;p + channels. Our results 
are A C p(rjK* Q ) = 0.17 ± 0.08 ± 0.01, A C p{vK* + ) = 
0.03 ± 0.10 ± 0.01, and A C p(riP + ) = -0-04i°j4 ± 0.01, 
all consistent with no asymmetry. 
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